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Abstract. Spherically symmetric brane black holes have tidal charge, which modifies 
. ^ ^ both weak and strong lensing characteristics. Even if lensing measurements are in 

^H' agreement with a Schwarzschild lens, the margin of error of the detecting instrument 

allows for a certain tidal charge. In this note we derive the respective constraint on 
&!)' the tidal charge of the supermassive black hole (SMBH) in the center of our galaxy, 

from the radius of the first relativistic Einstein ring, emerging in strong lensing. We 

find that even if general relativistic predictions are confirmed by high precision strong 
^ ' lensing measurements, SMBHs could have a much larger tidal charge, than the Sun or 

r — . neutron stars. 
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1. Introduction 

The Galactic Center is a highly dynamical region of the galaxy, however hard to resolve 
due to source confusion. Absorbing gas and dust is hiding the Galactic Center from 
Earth observers. Consequently, most of the knowledge about this part of our galaxy 
comes from observations in radio and infrared. A comprehensive introduction on the 
Galactic Center can be found in Ref . [1] . It has been estimated [2] that cold dark matter 
remnants of 1000 solar masses (M0) and compact star remnants of 1000 M© reside in 
the inner 0.01 parsec (pqll) region of the Galactic Center. The star population in the 
inner 0.04 pc, the SgrA* stellar cluster consists of B stars |3]. These are the remains 
from the dynamical process ejecting the hypervelocity stars of the Galactic halo [5]. In 
the distance range 0.04 — 0.5 pc Wolf Rayet and OB giant stars are found, all aged about 
6 million years [5]. 

The Galactic Center is dominated by a Supermassive Black Hole (SMBH), with 
mass A^B//=4.31 x 10^ M© and distance from the Sun tbh = 8330 pc, respectively 
[6]. The orbits of nearby stars depend on the mass and spin of the SMBH. These 
characteristics can be deduced from rigurous observations of stellar orbits [7],|8]. 

The first measurements of proper motions of stars within 2400 AU from the center 
of our Galaxy was published in Ref. [9j. A detailed review of the results of 16 years of 
monitoring stellar orbits around the SMBH using NIR techniques can be found in Ref. 
|10] . An animation of the orbits of individual stars around the SMBH has been produced 
from the images recorded between 1995 and 2011 by the UCLA Galactic Center Group 

Very recently a study of the stellar orbits around the Galactic Center by the 
NACO 9, SINFONI 10 and 11 observation programmes, performed with the Very Large 
Telescope (VLT), proved the existence of a gas cloud approaching the center of Sgr A* 
[6]. This can be seen from the redshifted H and He emission lines in the spectrum of 
the L' infrared band. The temperature of the gas cloud is 550 K, its luminosity 5 L© 
[6] and its mass of order 10^'* kg [12]. Table 1. of Ref. [6] summarizes the orbital 
elements of the cloud, implying an orbital period of 137 ±11 years. The velocity of the 
cloud increased from 1200 km/s in 2004 to 2350 km/s in 2011. The gas cloud began to 
disrupt in 2009 due to tidal shearing arising from the SMBH's gravitational force. Its 
shorter characteristic size is either marginally resolved or unresolved, smaller than 24 
milliarcsecond (mas). Its longer size is 46 ± 1 mas derived from the angular distance 
between the redshifted and blueshifted margins of the cloud. This is consistent with 
the photometrical width at half- maximum 42 ± 1 mas (2008) and 38 ± 1.6 mas (2011) of 
the cloud. It will reach the pericenter in 2013, at about 3140 times the event horizon of 
the SMBH. Only the two stars S2 and S14 have come closer to the black hole since the 
monitoring started in 1992 [6j. 

Light collected by all four telescopes of the VLT was successfully combined for 
the first time in 2011, using a new generation instrument in the VLT Interferometer, 

t 1 pc = 206260 astronomical units (AU) 
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PIONIER [13]. The design of another instrument, GRAVITY has reached completion 
as well [H] . This will be able to use the four telescopes as an interferometer in the near- 
infrared (NIR) band. The astrometrical accuracy strongly depends on the number of 
sources in the field of view, the source confusion being the main factor of the uncertainty 
of position measurement [15]. The statistical properties of the stars at the centre 
of the Galaxy make it plausible that at least three stars could be observable at any 
time in GRAVITY'S field of view [T3]. GRAVITY will perform astrometry with 12 
microarcsecond (/ias) precision in the K band up to 15 magnitudes [Hj. The margin of 
error of its measurements will be therefore 

A0 = 12 X 10"^ as . (1) 

From the numerous scientific applications of GRAVITY, described in Section 3 of 
Ref. [H], the observation of relativistic motions near the horizon of Sgr A* will be 
of uttermost importance. 

By such observations of NIR flares even the metric near the horizon can be 
determined [T6|. A submm-VLBI array should be able to actually resolve Sgr A* showing 
its event horizon as a shadow [17]. According to Ref. [7] the no-hair theorem can be 
tested. 

The planned accuracy of future radio instruments is compatible or even better (of 
order //as for the Square Kilometre Array [18] however for lensing phenomena we are 
primarely interested in NIR measurements. (Lensing observations in NIR are common, 
for an examle see Ref. [19].) 

With increasing measurement accuracy and the possibility of direct observations 
of the horizon of the SMBH at the Galactic Center the possibility of testing general 
relativistic predictions emerges. Even if a certain measurement confirms them, the 
measurement accuracy of a specific instrument will allow for certain margin of error for 
any additional parameter. From among the various possibilities we will test here the 
existence of a tidal charge, an imprint of the possible 5-dimensional nature of gravity in 
a brane- world scenario [20], [21]. A tidal charge emerges from the Weyl curvature of a 5- 
dimensional space-time in which our 4-dimensional observable world is embedded. Such 
a tidal charged black hole solution on a brane was found in Ref. [22]. The electric part 
of the Weyl curvature, non-standard model 5-dimensional fields, asymmetric embedding 
and a varying brane tension all could act as sources of the effective Einstein equation 
[23] . [21], [25]. In consequence unusual cosmologies may emerge [26], [27]. Remarkably, 
the Weyl curvature source term could replace dark matter [28], [29j . 

Light deflection and weak gravitational lensing by tidal charged brane black holes 
was investigated in Refs. [30], [31], [32]. One of the results of these investigations was 
a strict limit imposed on the tidal charge from Solar System measurements. 

In this paper we study strongly relativistic orbits in the Galactic Center. We focus 
on 1-loop photon trajectories, as depicted on Fig. 1. We summarize the basic equations 
concerning the one-loop null geodesies in spherically symmetric, static space-times in 
Section [2j The formation of the first relativistic Einstein ring, emerging due to strong 
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Figure 1. The trajectory of light along which the photons travel from the source S 
to the observer O, while turning around the lens L once, for q ~ Q, D]^ = 8600 pc, 
Dls = 10 pc. On panel (a) the whole geodesies is seen on logarithmic scale. The 
distorted region in the left and right side of the loop is magnified and shown on a 
linear scale on panel (b). 



lensing by the SMBH in the center of our galaxy is discussed in Section [3l under the 
assumption, that the metric is the tidal charged brane black hole |22]. Then we discuss 
the possibility that the radius of this Einstein ring varies in the margin of error of 
GRAVITY in the infrared astrometry regime and allow for a non- vanishing tidal charge 
in this range. 

We present our conclusions, which include the constraints on the tidal charge 
derived by this method, in Section |H 

2. Null geodesies and relativistic Einstein angles in spherically symmetric, 
static space-times 

The general spherically symmetric, static metric is 

ds"^ = gttdt^ + grrdr'^ + r'^dO'^ + r'^sin'^9 d^p^ . (2) 

Without reducing generality, we can restrict geodesic motions to the plane 9 = 7i/2. 
The (second order) radial geodesic equation can be replaced by the (first order) null 
condition 

where p is a parameter of the null curve. As the spherically symmetric, static metric 
does not depend on either of the coordinates t or ip, two constants of motion emerge by 
integration from the geodesic equations: 

r ^'^ (A\ 

i = 9..^ . (1) 

E ^ ft.| . (5) 

L is the specific angular momentum of the photon [33], its dimension being [length^], 
while Ec'' jG is the energy of the photon, the dimension of E being [length]. 



Black hole tidal charge constrained by strong gravitational lensing 



Then the rr term in Eq. (jH]) can be written in the form g^r {dr /dip) LP' / g"^^, while 
tt an 
becomes 



the tt and (pip terms as E'^/gu and L'^/g^^, respectively. Hence the radial equation 



E' drV L^ L^ 

9tt \dipj g^^ g^^ 

By reordering the terms, we obtain the equation characterising the trajectory r {ip): 

(7) 



■ /p2 „ \ni/2 

gifiifi j -'-' g(fiip , 



dr 

dip [ grr V -^^ -Qtt 

The sign differentiates between the incoming and outgoing parts of the path|§| The 
trajectory can be also expressed in terms of the dimensionless radial coordinate R = 
r/M. 

We will seek a null geodesic curve along which the photons travel from the source 
S to the observer O, while turning around the lens L once (1-loop orbit). The source, 
the lens and the observer are on the same coordinate line </? = 0, called the optical axis, 
more specifically L is the origin, S lies at (v? = 0, r = Dls), while O at {ip = 7i,r = Dl). 
Then the trajectory r {ip) of the photon is a decreasing function from r = Dls to 
some distance r = Tmin^ then increases to r = D^. This distance of minimal approach 
rmin = r (ipmin) IS a solutiou of the equation |3T] 

dr 

-r- {<^min) = . (8) 

dip 
For such a one-loop path the total change of the polar angle ip from 5 to O is tt + 27r. 
The function ip (p) is invertible either in the range n < ip < 27t (as there is only one r 
for each such ip), or in both monotonous parts of the path r (ip). (In these monotonic 
regions for any ip there is again only one corresponding r for each ip.) 

In what follows, we derive the angular radius 0^; of the first relativistic Einstein 
ring (the image created by all these 1-loop orbits). This is defined as the angle between 
the optical axis 5*0 and the tangent of the geodesic curve at the point O. Since at the 
point O the radial vector d/dR is parallel with the optical axis, we define the observable 
Einstein angle Qe as the angle between the vector d/dR and the tangent of the curve 
R (ip) at the point O. 

In the following derivation every vector field is evaluated at O. The scalar product 
(given by the metric ([2])) of the base vector d/dR and the tangent d/dp can be computed 
in two independent ways. First, from the general definition of the (Euclidian) scalar 
product of two vectors: 
d d 



d 

dR 


d 

dp 



dR dp " " 



^9rr{Dl)\I l^^y 9RRiDL)+ (^^y g,ADL)cosQE -(9) 

§ Eq. (8.5.6) of Ref. [SI] presents this equation in a reciprocal and integrated form, with the source 
and the observer placed at infinity, featuring r„iin instead of the impact parameter L/E; and employing 
the specific choice g^^ 



r 
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Second, from the decomposition of the tangent vector in the polar coordinate system 
we get: 

d d d fdR d du:> d \ dR d d dR 



dR d dip d , w^. ^ ^ v.^. , 
dR dp dR \ dp dR dp dip J dp dR dR dp 
Putting the right hand sides equal, we obtain 

Qe '■ = arccos 



.(10) 



-^\/gRR{DL) 



■g^j 9Rr{D 



arccos 



dR\ _ 
d^ J 9RR 



L) ^ 
-,1/2 



9w i^i 



dR\ 



:ii) 



S^j 9rr + 9w 

In the last step the parametrization of the curve R {p>) was eliminated by multiplying 
both the numerator and the denominator by dp/ dp (the angle between a vector and a 
curve must be independent from the parametrization). The right hand side of Eq. flTTl) 
is evaluated at i? = Dl/M (the location of the observer). 



3. The first relativistic Einstein ring and constraints on the tidal charge 



The tidal charged black hole [22] generalizes the Schwarzschild solution by allowing for 
a tidal charge parameter q in the metric function 



1 2M 
9tt = = -1 + 

9rT ^ 



q 



(12) 



of the spherically symmetric, static space-time (|2]). This emerges as a vacuum solution 
of the effective Einstein equation on a brane with q allowed to take either sign. For g < 
there is a horizon at M + (M^ — qY^'^, while for < g < M^ there are two horizons, at 
M ± (M^ — g)^/^. When q = M^ the two horizons coincide at M , while for q > M^ the 
metric describes a naked singularity. 

The trajectory of a photon, Eq. ([7]) in this space-time becomes 



dR _ 

dip 



R*-R' + 2R 



M2 



(13) 



Here the impact parameter L/E, normalized to M is dimensionless, similarly as R. The 
minimal dimensionless distance Rmin '■= ^mm/M (larger than the outer horizon of the 
lens) is a root of the polynomial emerging from Eq. ( 1T3l) by setting to zero the left hand 
side. 



The Einstein angle ( ITTl) then becomes 

me) 



e, 



arccos 



M 



2(^ 



M2 \dJ 



.(14) 



Next we assume a stellar source on the optical axis defined by the SMBH and the 
GALAXY detector, at a distance Dls varying in the range [10 pc, 100000 pc] and a 
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Table 1. Column 1.: the tidal charge, column 2.: the angular radius of the first 
relativistic Einstein ring, column 3.: the normalized radius of the horizon, column 4: 
the minimal distance Rmin = fmin/M, larger than the horizon in each case, column 
5.: the parameter L/ME, column 6.: the normalized tidal charge. 



q [W m2] 


Qe [l^-ds] 


th/M 


rrrdn/M 


L/ME 


q/M^ 


-1.815141 


38 


3.345 


AmA 


7.706 


-4.5 


-1.613458 


37 


3.236 


4.772 


7.496 


-4.0 


-1.411776 


36 


3.121 


4.612 


7.275 


-3.5 


-1.210094 


35 


3.000 


4.443 


7.043 


-3.0 


-1.008411 


34 


2.870 


4.265 


6.797 


-2.5 


-0.806729 


32 


2.732 


4.073 


6.535 


-2.0 


-0.605047 


31 


2.581 


3.866 


6.252 


-1.5 


-0.403364 


29 


2.414 


3.639 


5.943 


-1.0 


-0.201682 


28 


2.224 


3.385 


5.598 


-0.5 


0.000000 


26 


2.000 


3.090 


5.202 


0.0 


0.524374 


14 


n. s. 


1.446 


3.207 


1.3 



photon trajectory with one loop about the SMBH. Each value of L/ME determines a 
null geodesic curve by Eq. f lT3|) . however for most of the values the resulting curve does 
not reproduce the desired lensing geometry (e.g. the boundary conditions set by S and 
O). Therefore one has to 'fine-tune' L/ME to generate the specific curve which satisfies 
the conditions preset by the 1-loop orbit. We have done this in the following way. After 
setting the tidal charge q and the distance Dls, we choose some value for the L/ME (it 
turned out to be useful to chose it between 1 and 10). Then we evolve the differential 
equation ( TT3l) numerically with this value and calculate the change in the polar angle (p 
which emerges, while the photon travels from the distance D^s/M to Rmin then further 
to Di/M. For the choice L/ME = 10 this angular change was too large {A(p > Sir) for 
any q G [—5 M^, -1-1.3 M^] and Dls ^[10 pc, 100000 pc]. Therefore we have to reduce 
the value of L/ME, and repeat the procedure, until it yields Aip = 3tt (with 10^^^ rad 
accuracy). In this iterative way we obtain the dimensionless impact parameter for the 
considered 1-loop photon orbit. Finally we insert the correct value of L/ME in Eq. (TT^ 
and calculate the size of the first relativistic Einstein ring. 

The values of both L/ME and G^; found for Dls = 10 pc and q G 
[—5 M^, -|- 1.3 M^] are listed in Table [H The colums are (1): the discrete values 
of the tidal charge in the range studied, (2): the angular radius of the first relativistic 
Einstein ring, (3): the horizon radius normalized by mass (when applicable), or naked 
singularity (n. s.) for M^ < q, (4): the normalized minimal distance from the lens 
Rmin = ^min/M, (5): the dimeusiouless impact parameter L/ME, (6): the tidal charge 
normalized by M^. The g = line in the table reproduces the Enstein angle calculated 
for a Schwarzschild lens in Ref. [35], 0e(O, Dls) = 26 /las. 
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We have checked that by varying D^s the values of the Einstein angles change by 
less than 1/xas (this is how the agreement with the result of Ref. [35], calculated for 
Dls = 8600 pc, is acheved). Also note that the values of the tidal charge were varied 
such as the Einstein radius changes at most with the value ([1]), the designed margin of 
error of the measurements of the detecting instrument, GRAVITY, as compared to the 
Schwarzschild value. This allows for the limits q^ain and gmax given by 

QEiqn^a., Dls) = 6^,(0, Dls) - AG . (15) 

Thus a tidal charge falling in the range q G [—1.815, 0.524] x 10^° m^ will be consistent 
with measurements of the first relativistic Einstein ring generated by sources opposite 
to us with respect to the central SMBH. 

4. Concluding Remarks 

Bounds on the tidal charge of various astrophysical objects were derived earlier in the 
literature. For neutron stars a limit of |g| < 10^ m^ was establised in Ref. [36] from 
orbital models of high-frequency quasiperiodic oscillations observed in neutron star 
binary systems. From the constraint on the brane tension presented by Eq. (30) in 
Ref. [37] a weaker limit for negative tidal charges emerges in the following way. The 
junction condition (Eq. (35) in Ref. [37j) of the uniform density star and its exterior 
represented by the tidal charged metric associates a tidal charge to the limiting brane 
tension, provided the junction radius is known.. The latter is bounded from below by 
the compactness limit (the inequality (31) in Ref. [37]). This gives > g > —9.730 x 10^ 
m^. The strongest constraint for the Sun was found in Ref. [30] from the perihelion 
precession of the Earth, |g| < 6 x 10'^ m^. Light deflection measurements impose a 
milder restriction on the tidal charge of the Sun, |g| < 2.966 x 10^ m^ |31] . 

The question comes, whether experiments available in the near future targeting the 
observation of much larger objects in the Universe, where strong gravitational lensing 
could be relevant, would lead to other limits. Such tests concerning relativistic Einstein 
rings will be achievable by measurements of the high-precision designed instrument 
GRAVITY. We have studied the formation of the first relativistic Einstein ring in the 
tidal charged black hole geometry on the brane. For this we have specified the lensing 
geometry for the SMBH in the center of our galaxy, considered as a lens and for a light 
source opposite to us with respect to this SMBH, as a source. Following a similar logic 
to Refs [30], [31], thus assuming that strong lensing experiments will confirm the general 
relativistic predictions, but measurements error will still allow for some tidal charge, we 
have found a much larger value for the possible tidal charge of the SMBH in the Center 
of the Galaxy. Despite the instrument possibly confirming the Schwarzschild geometry 
within measurement accuracy, the range q G [—1.815, 0.524] x 10^° m^ of the tidal 
charge is still allowed, as the modifications induced by such a tidal charge fall within 
the designed instrument GRAVITY'S margin of measurement error 
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Table 2. Bounds on the tidal charge normahzed by mass square from observations in 
the Solar System (first column) [30) , constraints from orbital models of high-frequency 
quasiperiodic oscillations observed in neutron star binary systems [36 (second column) , 
constraints on brane tension and the compactness limit of neutron stars (third column) 
[57] , and finally, forthcoming strong lensing observations on the Galactic SMBH (fourth 
column), derived in this paper. Clearly, Solar System constraints are the strongest, 
neutron star bounds the weakest, while the neutron star binary and SMBH lensing 
constraints are of the same order of magnitude. 



object 


Solar System 


neuton star binary 


neuton star 


SMBH 


g/M^ ^^^ 


0.003 


2.339 


227.647 


4.485 



We have additionally checked that the second and third relativistic rings (with 
A(y9 = Stt, 77r, respectively) lead to similar results on the allowed range of g, as the first 
ring. We explain this by the relativitic rings being situated quite close to each other 
(within 0.03 yuas). 

Although the derived constraint on q is much weaker than those from neutron 
stars or Solar System observations, the dimensionless quantity q/w? is not very much 
different El In fact for this quantity the neutron stars constraints are the weakest, and 
the Solar System constraints the strongest, while both the neutron star binary systems 
and the SMBH strong lensing considerations presented in this paper gave comparable 
constraints. 

We conclude with the remark that even if general relativistic predictions on the 
Galactic SMBH are confirmed by high precision measurements, our investigations show 
that SMBHs could have a much larger tidal charge, than the Sun or neutron stars. 
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